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Abstract: 

We discuss the use of the variable V^miin which has been proposed in order to measure 
the hard scale of a multi parton final state event using inclusive quantities only, on a 
SUSY data sample for a 14TeV LHC. In its original version, where this variable was 
proposed on calorimeter level, the direct correlation to the hard scattering scale does 
not survive when effects from soft physics are taken into account. We here show that 
when using reconstructed objects instead of calorimeter energy and momenta as input, 
we manage to actually recover this correlation for the parameter point considered here. 
We furthermore discuss the effect of including W + jets and tt + jets background in 
our analysis and the use of V^min for the suppression of SM induced background in 
new physics searches. 



1 Introduction and motivation 



Since the startup and the following successful data taking of the LHC, the LHC exper- 
iments have already published a large number of result for exclusion limits for BSM 
physics, where for many BSM scenarios the actual limits within specific parameter re- 
gions have been strongly pushed to higher scales [1]. However, most of these analyses 
have been performed within specific models, and more generic variables, which provide 
information about the generic scale of new physics without additional assumptions 
about the decay topologies or decay chains, have not been fully exploited. Further- 
more, many variables which are currently proposed for mass or scale determination for 
new physics processes only make use of the transverse momentum of the event, thereby 
neglecting the information which can be obtained by additionally including the longi- 
tudinal information^. Examples for variables which make use of transverse momentum 
only are eg M5"2[4] or Mct [5]. On the other hand, some more traditional variables as 
eg invarant masses of composite objects as used in edges studies [6, 7, 8, 9] implicitely 
use all visible information, including the longitudinal momentum of the visible decay 
products. In [10] a new fully inclusive calorimeter-level variable, \/5min, was proposed 
which promises to give information about the hard scale of the underlying new physics 
processes without further assumptions or specification of the decay products, and addi- 
tionally makes use of the longitudinal information of the inclusive event. In that paper, 
the authors propose a conjecture which relates the peak position of \f§mia to the ac- 
tual rise of the hard new physics production cross section. However, this variable was 
subsequently shown to have a strong dependence on the soft physics in terms of ISR 
and underlying event. This is mainly caused by the fact that in this case, the energy of 
the additional soft particles equally enters in the calorimeter-level definition of \/5min, 
which then boosts the variable and its peak position to higher values with respect to 
the parton-level quantity. In [10], the authors tried to circumvent this problem by 
introducing a pseudorapidity cut in order to suppress the unwanted effects originating 
from soft physics from entering V^min- However, the introduction of the cut destroyed 
the correlation between peak position and hard cross section threshold which holds at 
parton level; this has explicitely been shown analytically for effects arising from initial 
state radiation (ISR) [11, 12]. In this case, the peak position is basically determined by 

the value of the pseudorapidity cut. Subsequently, y/Imin was promoted to V^min [13], 
using reconstructed objects at analysis level; for this variable, the correlation between 
its peak and the threshold of the hard production cross section was recovered, such 
that a determination of the hard scale of the BSM process was again made possible 
using experimentally accessible detector level objects. 

Apart from providing information about the hard scale of the underlying parton 
level BSM process, new variables can equally be used as cut parameters for SM back- 
ground suppression, and several of these variables have already made their way into 
the current BSM searches at the LHC experiments. It is therefore equally important to 
determine the use of \/Smin for SM background suppression. Although this constitutes 
a slightly weaker use of the variable per se, it is still an important issue to investigate, 

^Excellent reviews about different mass determination variables and their use, including advantages 
and disadvantages, have recently been published in [2, 3]. 
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especially as it has been proposed on a fully inclusive level and can therefore be applied 
without any further assumptions on the model or the specific decay chains and topology. 

In this report, we therefore investigate the properties of y/^mm at analysis level us- 
ing reconstructed objects. For this, we use a full sample for the mSugra point SPSla^ 
[15], which contains all strong production as well as all decay chains; for this parameter 
point, our sample therefore corresponds to the full data set which would be obtained 
from strongly interacting initial cascade particles in a realization of SPSla^. We include 
soft physics in terms of initial and final state radiation (ISR/ FSR), as well as a fast 
detector simulation. We test the correlation between the threshold of the hard process 
and the peak of V^min on parton level for the inclusive sample as well as exclusive 
dominant final states. In addition, we show that, for our sample, this relation can be 
regained on reconstruction level using quite simple analysis object definitions, and that 
major discrepancies between analysis and parton level quantities can be traced back 
to uncertainties in the reconstruction of tau jets. We equally comment on the power 
of the analysis level variable for SM background suppression, and compare to similar 
variables. All analyses are done for a LHC-like proton proton collider with a center of 
mass energy of 14 TeV and an integrated luminosity J C = 1 fb~ . 

The report is organized as follows: in Section 2, we briefly review the variable 
definition of V^min and define other kinematic quantities which were used in our study. 
In Section 3, we describe the data set we use in this study. Section 4 contains the 
comparison of parton and analysis level V^min, and Section 5 describes the inclusion 
of SM background and the use of V^min for SM background suppression, as well as a 
brief comparison with other (transverse) variables. We conclude in Section 6. 

2 Variable definition 

In this section, we will briefly review the original variable definition as well as its RECO 
level version; the interested reader is referred to [10, 13] for a more detailed discussion. 

In general, V^min is defined on an event by event basis as the minimal value for the 
partonic center-of-mass energy which is in agreement with the events' momentum 
configuration. It can be derived through a minimization process [10] as 

V^min (Minv) = \/ + + $1 + M^^ (l) 

^We are aware that this parameter point has recently been excluded by ATLAS measurements 
in the quark/ gluon plus missing transverse energy channel [14]. However, we here want to show 
that the parton level V^min can actually be recovered with sufficient accuracy from analysis level 
objects. Spectra which evade current exclusion limits typically exhibit higher initial cascade particle 
masses, and our arguments are generically not affected by the actual position of the particle production 
threshold. This is only important in the studies of background suppression presented in Section 5; here, 
a higher peak value for the BSM induced variable should actually even enhance the SM background 
suppression which can be obtained using Vimin- 

■^Gaugino-gaugino initial cascade states, which have not been considered here, would contribute an 
additional 5% to the total production cross section. 
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with 



VIS 



being the effective visible mass and the invisible total mass 

Minv = ^ m (2) 

invisible 

the sum over all masses of invisible particles. For completeness, we give the minimiza- 
tion procedure leading to Eqn. (1) in Appendix B. Note that Mvis and M^^^ are not 
defined equivalently, and especially that Mi^^ is not the Lorentz-invariant mass of the 
total invisible system, but rather the sum over all invisible particles rest masses^. In 
the definition of V^min, Mnv is therefore an external input parameter, as it is the only 
quantity which cannot be measured directly from experiment. Therefore, all results 
which are derived in the following sections have an implicit dependence on the value 
of Minv Throughout our study, we have usually set this to its "true" BSM value 
Minv = '2m^o{= 193.4 GeV). 

The translation to experimentally accessible quantities is then straightforward and 
gives 

. ^ - . I m - p2 A. J t'^ 

mm V mv 



(Mi„.) = ^E^ - P| + ^J$t, + Ml^. (3) 

Note that the use of transverse energy and momentum strongly depends on the defini- 
tion of the specific quantity; we define 

= -^T, $T = \^t\ (4) 
where E, P are the total energy and four momentum of all visible objects 

= Y^pl E = P\^r=(^j^\ (5) 

vis ^ ' 

and the z-direction defines the beam-line. In the original proposal, all visible quan- 
tities are taken from calorimeters, and soft background is suppressed by a cut in the 
pseudorapidity 77. Subsequently, for the correct value of Mjnv, true, a conjecture 

(Mnv, true) ~ V^th (6) 

J peak 

is empirically derived, which links the peak position of A/^min to the actual threshold 
^/s^Yi '^^ the hard matrix element process. However, in [11] it was subsequently pointed 
out that, using different values for the rj cut, the peak position for the calorimeter-based 
variable V^min could actually be arbitrarily shifted around; the same effect has been 
observed in [16], which applies the original calorimeter-based definition of V^min on 
our data sample. In answer to this criticism, new reconstruction and subsystem level 
variables were proposed in [13]. We here use V^min on an inclusive level using recon- 
structed objects, which basically corresponds to the RECO variable definition given in 



^For consistency, we here adopted the notation introduced in [10] for Minv and hope that the 
potentiahy misleading nomenclature does not cause confusion in the remainder of our discussion. 
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[13]. In our work, the suppression of effects from the parton shower has been achieved 
by quite simple object-level definitions given in Table 2. We will show that we obtain 
the parton-level v^min quite accurately in our sample, and that we indeed observe a 
similar peak of the RECO-level \f§m\a close to the production threshold for the correct 
input value of Mjnv In addition, we investigate the actual sources of discrepancies 
between the parton and reconstruction level -y/Smin in more detail. In our work, we 
equally present the first study of \/Smin as a variable for SM background suppression 
in BSM searches. 

In this study, we use the term "leptons" for all three SM lepton generations; in 
cases when we are concerned with tau leptons alone we will mention this explicitly. 
Equally, the tau jets at analysis/ reconstruction level are defined by the tau jet re- 
construction algorithm in Delphes [17] and differ from the parton level tau lepton by 
the four-momenta of the invisible tau decay products, specifically the associated third 
generation neutrino. In the following, we use the term "tau" for the parton level and 
"tau jet" for the reconstruction level quantity, which for an ideal reconstruction of 
the visible tau decay products four-momenta only differ by the four-momentum of the 
invisible decay products. 



3 Data sample and event generation 

In this report, we have made use of the BSM data samples which have been generated 
in the course of the 2009 BSM Les Houches mass determination study; first results 
using these data for studies of various mass determination methods were presented in 
[16]. We use a SUSY spectrum for the point SPSla, where the spectrum was generated 
with the spectrum generator SOFTSUSY [18]. Parton level events have been generated 
using Madgraph [19, 20] for the generation of the heavy initial cascade particles (i.e. 
the squark-squark, squark-gluino, and gluino-gluino initial states). The heavy pair- 
produced particles have then been fully decayed according to the respective branching 
ratios into all possible decay products using Bridge [21] within the Madgraph frame- 
work; we therefore consider a complete sample for this parameter point, which contains 
all possible final states. The SM background has been generated using Alpgen [22]. 
For parton shower and hadronization, we used Pythia [23], where the parton shower 
evolution follows the Pythia 6.4 default, ie is ordered with additional modifications 
to guarantee color coherence, as well as matrix element corrections where these are 
available (Pythia switches are given in Appendix C). The detector simulation has 
been performed with Delphes [17] in its default mode. For specific input parameters 
and setups, we refer to the specifications which can be found in the data base for our 
samples [24]. Data analysis as well has fitting has been done within the ROOT [25, 26] 
framework. All our results have been obtained with a data sample for a center-of- 
mass energy A/S'hadr = 14TeV and an integrated luminosity j C = lfb~^. Table 1 
lists the production cross sections for the hard 2 — !■ 2 process; these numbers were 
obtained using the Madgraph parton level 2 — )■ 2 production cross sections, with the 
electroweak scale spectrum obtained from SOFTSUSY, convoluted with PDFs to ac- 
count for the parton to hadron transition for the incoming states. For this study, we 
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XIX2 


2^2 


qq 


6.56 


Q9 


19.96 


99 


4.53 



Table 1: SPSla production cross sections in pb for pp — X1X2 using Madgraph 
2 — >■ 2 parton level production cross sections, convoluted with PDFs, for a hadronic 
center-of-mass energy of 14TeV. CTEQ6L1 PDFs [27] were used. 

restrict ourselves to the leading order predictions for the hard process in both signal 
and background simulation^. 

4 Parton and analysis level V^min 

As already discussed in Section 2, the variable V^min has undergone several develop- 
ments since its original proposal. Initially defined as a calorimeter-based variable, it 
was shown to be quite sensitive to effects of soft physics for the respective processes. 
Especially the original merit of this variable, namely the correlation of the peak position 
and the threshold of the heavy pair-produced particles at the beginning of the decay 
chain, is strongly influenced by the soft physics of the event. The original suggestion 
of the authors was to introduce a cut on the pseudorapidity; however, the authors in 
[11, 12] have shown analytically that the position of the peak position in this case is 
completely cut- value dependent; similar results have been observed in [16]. 

In this study, we show that, if y/^mm is defined at analysis object level rather than 
on calorimeter level, the parton level variable V^min can be reconstructed quite well 
using simple object definitions. This recovery of the parton level peak position using 
a reco-level variable for both inclusive and exclusive final states, however for different 
parameter points, have equally been presented in [13]. For our sample, we equally 
observe that the conjectured correlation between the rise of the hard scattering event 
cross section and the peak position of y/Imm holds; however, we want to emphasize that 
this is on the level of a conjecture which has not been systematically studied or proven 
on an analytic level, although some preliminary studies indicate a kinematic origin 
which emerges after the convolution with PDFs^ [28] . Therefore, even on parton level, 
it is currently unclear whether this conjecture necessarily holds for all BSM parameter 
points and scenarios. Equally, this conjecture only holds for a correct input value of 



fully differential study including NLO contributions to account for cut effects would require 
2 — )• n event generators for both BSM signal and SM background, which additionally include the 
matching of parton shower and NLO contribution; although fast progress has been made in this field 
for SM processes, no fully differential higher order BSM generator is currently publicly available. 

^This result has been obtained with a unit matrix element as well as unit PDFs; in this case, the 
peak position of v^min arises from the lower PDF integration boundary following from the kinematic 
lower limit which guarantees that y^part ^ V^throshoW- More realistic scenarios with non- uniform 
PDFs and matrix elements are currently under investigation. 
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object 


Delphes predeRnition 


additional requirement 


electron/ position 


\ri\ < 2.5 in tracker, px > lOGeV 


isolated 


muon 


\7]\ < 2.4 in tracker, px > lOGeV 


isolated 


lepton isolation criteria 


no track with px > 2 GeV 


no track with px > Q GeV 




in a cone with dR = 0.5 


in a cone with dR = 0.5 




around the considered lepton 


around the considered lepton 


n leptons 





exactly n isolated leptons 






at detector level 


taujet 


PT > 10 GeV 




jet 


PT > 20 GeV 


PT,jet > 50 GeV, Iryljet < 3 




CDF jet cluster algorithm [29], 






R = 0.7 




Missing transverse 




> 100 GeV 


energy 







Table 2: Physical object definitions in terms of the single objects pseudorapidity rj, 
absolute value of transverse momentum px, and (jet) cone radius R for analysis level 
objects on detector level. We basically adapt the Delphes predefinitions, with slightly 
more stringent requirements for isolated leptons and jet definitions. We equally set a 
lower limit i?™'^^ > 100 GeV for events with missing transverse energy. 



In the following, we will compare quantities derived on the parton level with the 
same quantities which have been derived from analysis level objects. For the identifica- 
tion of the former, we consider the hard process, i.e. our data sample after the complete 
decay to SM particles and the LSP, but before the parton shower, hadronization, and 
detector simulation. All particles are considered as visible apart from neutrinos and 
the LSP. The invisible total four-momentum is then the sum of the latter particles' 
four-vectors 

pparton _ 

-"^ invis / J P^y 

and the same holds for the missing transverse momentum. At analysis level, we re- 
quire all physical objects to fulfill the object definition requirements given in Table 2 
on detector level; these object definitions basically follow the Delphes predefinitions, 
where we introduced slightly more stringent requirements for lepton isolation and jet 
criteria and equally set a lower limit of 100 GeV on the total missing energy^. Visi- 
ble and invisible quantities are then defined according to Eqns. (4) and (5) in Section 2. 

We first study the variable V^min for a complete inclusive sample, i.e. we sum over 
all final states of the hard process. Our main results are shown in Figure 1, where 
we compare the true parton level v^min, reconstruction level A/^min as well as the 
original calorimeter based variable with and without a cut in pseudo rapidity rj, as 
originally suggested in [10]. We see that the parton level y/Imin peaks quite close to 
the actual heavy particle production threshold as suggested in [10]; equally, we observe 

^These cuts closely follow cuts used in the SUSY analysis studies in [8, 30]. Due to the relatively 
high PT jet cuts, together with a high E™^^^ cut and lepton isolation criteria, we expect minimum bias 
events to be sufficiently suppressed ([31], as well as section 6.1 in [32]). 
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that the same variable from reconstructed objects again peaks close to the threshold, 
but is shifted to slightly lower values with respect to the parton level quantity. We will 
comment on this in more detail below. In contrast, the pure calorimeter based variable 
exhibits a peak position at quite high values and can therefore not be used for a scale 
measurement of the new physics process. Restricting the contributions to calorimeter 
energy deposits with a minimal pseudorapidity improves this behavior and brings the 
peak closer to lower values; however, this approach suffers from the drawbacks pointed 
out in [11, 12]. 

For a more accurate determination of the peak position and a viable assessment 
of the error in its position, we fit the V^min distribution with a Gaussian around its 
peak, where we use the largest fit region which is still in agreement with x^/d.o.f ~ 
C(l). Specifically, we use 600 GeV < v^^in < 1400 GeV and 400 GeV < v^n,in < 
1400 GeV to determine the parton level and analysis level peak positions respectively. 
We then obtain 

parton level V s^^^ : (1152 ± 4) GeV 
analysis level V^min : (1083 ± 4) GeV 

We see that the reconstruction level variable for the overall sample peaks close to 
the "true" maximum of the parton level variable, the difference being (9(100 GeV). 
In order to pin down the major sources of this shift, we have performed detailed 
studies for specific final state signatures; we will discuss this in more detail in Section 
4.1. To summarise the result of this section, we observe that, in our sample, larger 
shifts in the peak positions stem from processes with one or more leptons in the final 
state. One source of this is the imperfect reconstruction of tau jets from parton to 
analysis level objects. We can test this by taking an "idealistic" approach, where 
we use the parton level four-vector values for taus in the analysis level objects; this 
simple "gedankenexperiment" trick, where we assume a perfect reconstruction of tau 
jets at analysis level, significantly reduces this difference, cf. Figure 2. While such a 
requirement is in fact not possible in reality, it however shows that our (quite loose) 
lepton definitions and resulting poor tau reconstruction are a major source of this 
shift, and more dedicated algorithms might further reduce this discrepancy. Fitting 
the "new" analysis level distribution within the range 600 GeV < v^min < 1400 GeV, 
we obtain 

analysis level V Sj^^^^ , r = Tparton : (1163 ± 4) GeV 

and we see that the discrepancy with the parton level value of (1152 ± 4) GeV reduces 
to the permill level, cf. Fig. 2. The average heavy particle threshold in our sample is 

true (average) (mi + 1712) ■ 1146 GeV 
which again agrees with the parton level value of V^min on permill level within the 
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Figure 1: Sum of qq, qg, , and gg initial states, where qg initial states dominate. All 
final states which fulfill object definitions from Table 2 are included. True y/S (red; 
solid), parton level VSmm (blue; dashed), analysis level VSmin (green; dotted), y/s^^^ 
using calorimeters (pink; dash-dotted) and same with an |?7| < 1.4 cut (black; dash- 
dot-dot-dotted). V^hadr = 14TeV, J jC = lfb~^; corresponds to 31050 events. Shift 
between parton level and analysis level peak is about 70 GeV. The calorimeter based 
distribution without a pseudorapidity cut exhibits a peak at much larger \/Imm values. 

error bars ®. In addition, the object definitions in Table 2 equally allow for an adequate 
recovery of the parton level distribution shape, and, more specifically, we are able to 
suppress distribution tails for higher -\/Smin values appearing in the reco-level definition 
of this variable in [13]. A breakdown in terms of pairs of initially produced particles 
prior to the cascade decays is given in Table 3. 

For illustration purposes and completeness, we also investigate the parton level 
V^min dependence on the input value for Mi^v', similar results have already been pre- 
sented in [10] and [16] for the analysis level quantity. For R-parity conserving SUSY 
scenarios, as considered in this study, this corresponds to the guess of the LSP mass, 
as in this case Mi^v = 2mLSP when neutrino masses are neglected. Figure 3 shows 
the shift of the parton level V^min distribution for different input values Mjny. We see 
that a variation of the M[^^ mass leads to a shift in the peak of a similar magnitude. 
Therefore, we again emphasize that the results presented in this study concerning the 
correlation of the peak of the V^min distribution and the hard scale of the underly- 
ing production process have indeed an implicit dependence on the correctness of the 
guessed input value for Mjnv, as already discussed in the original proposal of this vari- 
able [10], and therefore generically only allow for a measurement of the hard scale as a 
function of this variable^. As before, the parton level distributions and peak positions 

^In this work, we only want to demonstrate that the approximate peak position of the parton level 
variable can actually be obtained from analysis level objects; for more dedicated analyses, the peak 
position could also be determined by other means, eg. a fit to a more variable-specific function. 

^Several other widely used variables, as eg the original definition of Mt2 [4], equally exhibit a 
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Figure 2: As Fig. 1, where the hard matrix element tau four-vectors were used for the 
analysis level V^min- With the differences due to tau identification at the analysis level 
removed, parton and analysis level peak positions agree within error bars. Explicit 
numbers are given in Section 4. 



initial cascade particles 


threshold 


r^peak;parton 


/^peak;ana 


/^peak;ana,r = Tp 




[GeV] 


[GeV] 


[GcV] 


[GeV] 


9Q 


1163 


1195 ± 5 


1101 ± 6 


1204 ± 6 






1170 


1035 


1170 


qq 


1046 


1042 ± 5 


1012 ± 8 


1088 ± 8 






1080 


1035 


1080 


99 


1215 


1257 ± 7 


1150 ± 7 


1241 ± 6 






1260 


1170 


1260 



Table 3: Peak positions for separate heavy initial cascade particles for parton level 
(Vlj^'if '''"'°") and analysis level (v/I^lf quantities as well as analysis level quantity 
for idealized tau jets (VifX''''"''''"""'') i^i GeV, where the value in the respective first line 
arises from a Gaussian fit around the peak, while the second corresponds to the more 
simplified definition of the peak position by maximal number of bin entries. In addition, 
we give the average threshold value y/s^^ = (mi + 777.2) for each sample, where mi 2 
are the masses of the heavy initial cascade particles. We see that the peak position 
from both peak position definitions are close to the actual thresholds; in addition, the 
effect of imperfect tau reconstruction account for an approximate shift (9(100 GeV) for 
all initial state pairings. 
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Figure 3: Parton level Ai/Ipart (red; solid), and parton level V^min depen- 
dence on the input value Minv for varying values. Shown are results for 
Minv = GeV (dark green; long dashed), 200 GeV (black; solid), 400 GeV (pink; dotted), 
600 GeV (green; dash-dotted), 800 GeV (dark blue; dash-dot-dot-dotted), 

1000 GeV (light blue; short dashed). The corresponding peak positions for increas- 
ing Minv input values as given above, using a maximal bin definition, are obtained as 
^peak; part ^ ^^^^^^ ^^^^^ ^^^^^ ^^^^^ ^^^^^ ^^^^^ before, the analysis level 

peak positions using the parton level tau leptons coincide with the parton level peaks 
(not shown here). 

could be reproduced using analysis level objects in the idealized version, ie replacing 
the analysis level tau-jets with parton-level tau leptons in V^min- 

4.1 Signal based searches 

In this section, we consider the variable V^min for several exclusive final states. The 
(parton level) dominant decay modes of our sample are given in Table 4. 



Most of the dominant final states can be tracked down to a couple of competing pro- 
cesses, and can be broken down to the following parton-level decay chains 

• g^, 3 jet channel 



qng qmrnRXi QBqRQRXiXi (90%) 

dependence on an input value for the LSP mass. 
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final states, hard matrix element 


main source 


* nara 


-^ana 


leptons, 3 jets 


Q9 


4480 


14247 


2 leptons, 3 jets 


qg{97%) 


4020 


2092 


1 lepton, 3 jets 


qg{99.99%) 


3740 


5282 


2 leptons, 2 jets 


qq 


1776 


2745 


1 lepton, 2 jets 


qq 


1366 


6997 



Table 4: Number of events for dominant parton level decay modes, characterized by 
specific visible final states, on parton level (A^hard) and at analysis level (A'^ana)- At 
analysis level, the jet number requirement for event selection is changed from an exact 
equality to a minimal number of jets. If not stated otherwise, the main source provides 
all events with a specific signature on the parton level. Examples for dominant decay 
chains leading to the specific parton-level final states are given in Section 4.1. 

• g^, 3 jet 2 lepton channel 

qL~9 -> qLT+T- xWrQrXi (27%) 
qR~9 -> qRX>bT+r-xl (22%) 
qR~9 ^ qRxU'LQ'Lr+r-xl (17%) 

~ ~ , / ~0 // -/// ~0 / 1 ^7a/ \ 

qL~9 ^ q'LTi^rXibtTU^Xi (17%) 



• g^, 3 jet 1 lepton channel 

qRg -> qRxlbirUrxl (45%, BR ~ 0.09) 

qR~9 qRXiQLq'LTi^rXi (30%, BR ~ 0.06) 

qL~9 ^ q'Lru^xWkQRXi (25%, BR ~ 0.05) 



• gg, 2 jet 2 lepton channel 

~ -/ ,11 ~0 III ~0 tnr(y/\ 

Ql qL Ql'^^tXi Ql '^^t Xi (36%) 

Qr^'l ^ qRxU'Lr'-r'^^ (64%) 



• gg, 2 jet 1 lepton channel 

QR^L ^ QRxULri^rXl (100%) 
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At the analysis/ reconstruction level, we here require to have a minimal jet multiplic- 
ity, which leads to much larger event numbers especially for signatures with a smaller 
number of leptons. We equally do not apply any dedicated additional channel-based 
cuts. 



Figures 4 and 5 show the true a/S, parton level V^min, reconstruction level V^min as 
well as the original calorimeter based variable with and without a cut in the magnitude 
of the pseudo rapidity |?7| < 1.4 for several explicit final states. We observe a similar 
behavior as in the overall sample, cf. Fig. 1: the parton level variable peaks around the 
actual production threshold, while there is a shift to lower peak values for the analysis 
level quantity. In order to understand the origin of this shift, we investigate this for 
a final state which initially exhibits a large difference between these quantities. We 
consider the 2 jet 1 tau-lepton channel, where originally the y/Smm peak positions differ 
by about 170 GeV. From Eq. (3), we see that the definition of V^mm depends on the 
following independently measured quantities: 

-^vis) Pz, I^tI = Pt = Et = t'T- 

In order to investigate the origin of the shift between the parton level and analysis level 
peak positions, we therefore consider each of these variables separately and plot the 
difference between the respective parton level and analysis level quantity; the results 
are shown in Figures 6 and 7. We observe that, while the differences between parton 
and analysis level Pz, Pt basically peak around zero, there is an average discrepancy 
~ 50 — 100 GeV between the parton and analysis level total visible energy. However, 
this discrepancy is accounted for by the fact that when changing from parton to hadron 
level, we replace the (visible) parton level tau by the (visible) tau-jet and the (invisible) 
tau-neutrino: 

In this transition, we equally shift the four momenta of the tau neutrinos from the 
visible to the invisible contribution of the definition of A/Imin (cf. Eq. (3)): 

pP^'' = ... + v^' + MFJ^'' = 



VIS 



-'vis ■■■ ' -f^^jet '■■■■> ^"mv ^"inv ' 



m,. 



We consider the neutrinos to be massless; therefore, we can leave the sum of all invisible 
particles' masses Minv unchanged. As the original variable definition of V^min and the 
subsequent correlation in Eqn. (6) only depend on the heavy initial cascade particles, 
but not on the actual number of visible and invisible decay products, the observed 
change in the visible energy due to the escaping neutrinos at analysis level should then 
be compensated by associated changes in Pz-, Pt on an event by event basis, leading to 
a similar peak behavior of the distribution at parton and analysis level. Here, in 

order to assess the overall impact of this shift and a possible poor reconstruction of the 
tau decay products^'', we perform a gedankenexperiment and change into a more ideal 
world where we idealistically reverse the analysis level tau-jet reconstruction and take 



-'^^Note that, in case the shift cannot be explained by poor reconstruction of the decay products alone, 
this equally opens the window to a possible topoiogy-dcpendcncc of V^min! we thank K.Sakurai for 
pointing this out. 
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Figure 4: as Fig. 1 j£ = Ifb ^: Left figure for exactly leptons in the final state, 
3 hardest jets, 14249 events), and shift between parton level and analysis level V^min 
is about 100 GeV. Right figure for exactly 2 leptons in the final state, 2 hardest jets 
(2745 events). Here, the shift between parton level and analysis level y/^mm is about 
200 GeV (reduces to 100 GeV if parton level tau vectors are used) 



final state 


/^peak;parton 

vs^in [GeV] 


Vs^in [GeV] 


^pcak;ana,r = Tp 
VSmin [GeV] 


0/3j 


1190 ± 5 


1072 ± 6 


1072 ± 6 


2Z3j 


1271 ± 8 


1128 ± 8 


1257 ± 8 


lr3j 


1204 ± 7 


1123 ± 8 


1210 ± 8 


2/2j 


1231 ± 7 


1001 ± 7 


1105 ± 6 


lr2j 


1157 ± 7 


990 ± 7 


1031 ± 8 



Table 5: Comparison of peak positions from Gaussian fits for parton level V^min 
{VsZf''""'"'') and analysis level V^min (v^mif'""") for specific final states, specified by 
the number of visible final state leptons (/), jets (j), and r-leptons (r), corresponding 
to dominant decay chains in the complete SPSla sample. Values for the peak position 
of the analysis level quantity with perfect tau jet reconstruction (v^lJ"^^""'""'''^ " ^'') are 
also given. For most final states, the effect of the peak shift due to imperfect tau jet 
reconstruction is (9(100 GeV). 



the parton level tau four-vectors for the analysis level variable. In this case, the shift 
in the peak position of y/lmm reduces to roughly 130 GeV. An alternative though less 
sophisticated way to determine the peak position is to consider the bin which contains 
a maximal number of entries; using this definition of the distribution peak position, 
the original shift between parton level and analysis level V^min reduces from 200 GeV 
to 90 GeV if the parton level tau vectors are used at analysis level. A similar study 
for the 2 tau lepton 2 jet channel, which originally equally exhibits a quite large shift 
between the peak positions, shows that the effect of tau misidentification is O{100GeV) 
for both peak position definitions, reducing to ~ 100 GeV in both cases when tau 
misidentification is removed. A similar effect can be observed for other specific final 
state signatures, cf. Table 5. 



Although we still obtain a quite large shift for specific final state signatures, we have 
seen that, when using parton level taus for the analysis level observable and therefore 
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sqrts(min) (GeV) sqrts(min) (GeV) 

Figure 5: as Fig. 1 J C = lfb~^, exactly 1 tau lepton in the final state, 2 hardest 
jets (3260 events) Shift between parton level and analysis level peak is 170 GeV(left) 
and reduces to 130 GeV(right) when parton level tau vectors are used for analysis level 
objects. 




Pt(parton)-PI(ana) (GeV) Elau(parton)-Etau(ana) (GeV) 



Figure 6: Difference between parton level and analysis object level total transverse 
momentum (left) and tau jet energy (right) for the 1 tau 2 jet channel. While the 
shift between the two values for the transverse momentum peaks around around zero, 
the shift between the parton level tau and analysis level tau jet energy is quite large, 
due to the escaping neutrino in the tau jet reconstruction. The difference in the Pz 
distribution (not shown here) exhibits a similar behaviour as the Pt distribution. 
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Figure 7: Difference between parton level and analysis object level total visible energy 
for the 1 tau 2 jet channel. Left side shows the real difference, while on the right hand 
side analysis level tau jets were replaced by parton level taus. While we originally 
observe a large shift between the two values, originating from the escaping tau neutri- 
nos on reconstruction object level and with a peak on the order of (9(200 GeV), the 
distribution of the difference peaks around zero when parton level tau four- vectors are 
used for the calculation of the analysis level observable. 



suppressing possible effects from poor tau reconstruction, the inclusive sample peaks 
at the same value for both parton and analysis level v^mm distributions, cf. Fig. 2. 
We therefore conclude that, with correctly identified analysis level objects, the peak 
position of the parton level y/Imm can indeed be reconstructed from generator level 
measurements^^; however, we want to emphasize that the correlation between the peak 
position and the actual heavy particle production threshold only exists in the form 
of a conjecture which lacks a rigorous proof. In case the conjecture proves to hold 
in all cases, the analysis level y/^tnm variable indeed gives a quite easy grasp on the 
threshold of the new physics pair-produced particles. Although this analysis was done 
in a specific scenario, where only certain initial heavy particle spin states are allowed, 
we saw that our conclusions hold for all possible spin combinations we considered. As 
our study relies on purely kinematic variables, we are therefore confident that these 
also hold for other spin combinations both for the heavy initial pair-produced particles 
as well as the particles in the decay chains, i.e. especially for other (also non-SUSY) 
BSM scenarios. 



4.2 Comment on additional soft physics effects 

The data set used in this study contains soft physics in the form of initial and final state 
radiation as described in Section 3, but no simulation of underlying event or pileup. 
However, the criticism which was expressed by the authors of [11] exactly concerns the 
dependence of soft physics in terms of ISR, which has been addressed in this work. 

^^We want to point out that the reconstruction level objects in Table 2, through their definition 
by pt and 77 cuts, still depend on these two parameters; therefore, a recovery of the hard scale from 
reconstruction level objects will always be obstructed by an implicit dependence on the cut values 
in the analysis object definitions. However, in contract to the calorimeter-based variable and the 
cut in pseudorapidity originally proposed in [10], we here use object level definitions which are more 
optimized to the reconstruct the hard scattering event. We thank B. Webber for bringing this point 
to our attention. 
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Additionally, soft physics can enter in the form of minimum bias events, underlying 
event and pileup. We believe that the cuts in Table 2 are sufficiently hard enough to 
suppress minimum bias events ([31], as well as section 6.1 in [32])). Underlying event as 
well as pileup effects can still distort the overall result for the peak position; however, we 
believe that these issues should be pursued in an experimental study, in combination 
with a collaboration internal full detector simulation. We can give a first estimate 
of the effect of underlying event fake Pt contributions by adding AP|?''*° = 10 GeV, 
which corresponds to a conservative upper limit of the average Pt from the underlying 
event [32, 33], in the definition of in Eq. (3); in this case, the best fit value from 
400 GeV < A/Jmin < 1400 GeV is again given by 

\fs2l^{/\P^^^'' = 10 GeV) : (1082 ± 4) GeV 

which completely agrees with the value without the addition of APl?*^*^. Changing the 
fake additional transverse momentum to AP^^*^ = 100 GeV leads to the result 

Vs^^^iAP^^^^ = 100 GeV) : (1093 ± 4) GeV. 

In fact, as M-^v and appear in the same form in the definition of y/^mm, the generic 
effects of additional fake P^s can be estimated from Figure 3. From underlying events, 
^pftikc < j^gg QgY^ g^^^j therefore estimate the uncertainty related to underlying 
event fake transverse momentum to be generically much smaller than the tau recon- 
struction effects discussed above. A more realistic investigation of these experimentally 
dominated effects, which should include a full detector simulation, is beyond the scope 
of this work^^ . 



5 SM background 

In this section, we investigate A/Imin when SM background is included, as well as its 
use for the reduction of SM background in new physics searches. As an example, we 
consider W + jets and ti + jets background. Due to the large cross sections, we applied 
an additional filter in the generation of the SM data sample, cf. Table 6. The cross 
sections after these additional cuts are given in Table 7. As before, we investigate the 
peak position for A/Smin at parton and analysis level when using the true BSM value 
Minv = 2 m^o. From Fig. 8 we see that in the total number of events after the cuts the 
peak structure disappears when no further SM cuts are applied. However, assuming 
an accurate enough (data or Monte Carlo driven) background subtraction, the peak 
structure is clearly visible again and much larger than the statistical error, cf. Fig. 8. 
A similar behavior is observed when we vary the input variable Mi-^v'. Fig. 9 shows the 
behavior for Mjnv = GeV, 1000 GeV respectively after SM background subtraction; 
we see we obtain a clear BSM signal. We therefore conclude that, at least for the 

-'^^Simulations using more recent tunes for a 14TcV LHC, as eg the Pythia C4 tunc [34], point to 
additional AP^^° ~ 30 — 40GcV for similar/ less stringent object definition values [35]; however, 
following the above considerations the induced error then is certainly ^ GeV, which again corresponds 

to a relatively small uncertainty in the determination of V^min ■ 

^■^In fact the experimental collaborations are already applying algorithms to subtract Et due to 
underlying event; cf eg [36]. We thank S. Wahrmund for bringing this to our attention. 
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''^Icptons 


\f T,mm 


< 2 


80GeV 


2 


40 GeV 


> 2 


OGeV 



Table 6: Additional filters on magnitude of the total missing transverse momenum \fT\ 
applied for SM background generation, depending on the number of final state leptons 
'^leptons- Leptons are required to obey the cut criterium \pt\ > 5 GeV for the magnitude 
of the transverse momentum and |?7| < 3.2 for the magnitude of pseudorapidity. 



'^jcts 


W + n jets 


tt + ra jets 







75.2 


1 




48.5 


2 


188 


20.0 


3 


53.6 


6.3 


4 


12.6 


1.4 



Table 7: Cross sections in pb for SM background with Alpgen; filters in Tab. 6 were 
applied in the generation stage. 




500 1000 1500 2000 2500 3000 3500 4000 4500 500 1000 1500 2000 2500 3000 3500 4000 4500 

sqrts(min) (GeV) sqrt5(min) (GeV) 



Figure 8: Analysis level Vlmin after a cut VImm > 700 GeV. Dominant six SM back- 
grounds after cut {W + 2j, W + 3j,W + 4 j, tt, tt + Ij, tt + 2 j) are included. Left: 
SM+BSM (green, dotted; 136834 events) and SM only (black, solid; 108017 events). In 
the sum and without further suppression cuts, the peak structure disappears. Right: 
Difference between (SM+BSM) and (SM). Assuming the SM background is well-known, 
the peak structure of the BSM signal is recovered. The difference is much larger than 
the statistical error. 
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Figure 9: Difference between (BSM+SM) and (SM) for Mjnv = GeV and a cut 
V^min > 500 GeV (left; 29815 events) as well as M-mv = 1000 GeV and a cut V^min > 
1500 GeV (right; 27802 events). Assuming the SM background is well-known, the 
peak structure of the BSM signal is recovered. The difference is much larger than the 
statistical error. 



parameter point studied here, with SM background being well-known, y/Imm can be 
used as a BSM discovery variable and that for a true input value of Mi^v, the scale 
of the new physics can be derived from the peak of both parton level and (properly 
defined) analysis level quantities. 



We furthermore assess the use of V^min as a cut variable for SM background suppres- 
sion. For this, we investigate the position of the peak for the different SM background 
channels considered here, where we again use Mjnv = 2m^o. The respective values are 
given in Table 8. For a first estimate of these positions, we do not need to perform 
a more sophisticated fit, and we therefore follow the simplified approach by defining 
the peak positions according to the bin which has the maximal number of entries. We 
see that the most dominant SM background channels have distribution peaks around 
500 — 700 GeV, while the BSM signals peak at higher values. We therefore apply two 
different cuts of V^min > 700 GeV and V^min > 800 GeV on all samples; the cross 
sections after these cuts are summarized in Table 9. We see that, for both cut values, 
while we only cut out around 10% of the BSM signal, the dominant SM channels are 
suppressed by factor 3-6. We therefore conclude that V^min can easily be used as a 
variable for SM background suppression, even for wrong guesses for the total invisible 
mass Minv. In the previous sections, we discussed how in our sample the peak of the 
V^min variable is correlated with the real threshold for the hard production cross sec- 
tions only if the correct value input for Mjnv is used, cf. Eqn. (6). For the background 
suppression, however, a correct guess or estimate of this value from other sources is 
not necessary, and we equally obtain a good SM background suppression with wrong 
input values for Min^^^. 



^^We want to remind the reader that the same value of Mjnv for the calculation of v^min needs to 
be used in both SM and BSM samples; the correlation with the threshold however only holds for the 
sample with the equivalent correct Minv We thank K. Matchev for reemphasizing this point. 
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/^maxbin;parton 


/^maxbin;analysis 


0", Vimin>700GeV 


0", Vimin>800GeV 


process 


a 


V S ■ 
* iiiin 


V S 

^ mill 




[pb] 


[GeV] 


[GeV] 

L J 


[pb] 


[pb] 

Li- J 


W + 2j 


188 


450 


405 


29.03 


19.07 


W + 3j 


53.6 


630 


585 


23.47 


17.51 


W + Aj 

r J 1 J. J 


12.6 


900 


765 


9.03 


7.63 


y W + iets 


254.2 






61.53 


44.21 


ti 


75.2 


540 


450 


11.50 


3.35 


ft A- I i 


48.5 


675 


585 


20.71 


14.41 


tt + 2j 


20.0 


900 


720 


14.28 


11.68 


tt + 3j 


6.3 


1215 


900 


5.31 


4.85 


tt + Aj 


1.4 


1530 


1215 


1.29 


1.24 


E t + jets 


151.4 






53.09 


35.53 


qq 


6.56 


1080 


1035 


5.59 


5.03 


99 


4.53 


1260 


1170 


4.38 


4.20 


m 


19.96 


1170 


1035 


18.85 


17.82 


BSM 


31.05 






28.82 


27.05 



Table 8: Cross sections a for SM background processes and "\/Smin maximal bin po- 
sitions for parton level (Vi:::r,f and analysis level (Vll^ff with standard 
analysis object definitions only; Min^ = 2m^o. Last two columns give cross sections a 

after \/Smin cuts respectively. After a minimal analysis level cut on V^min, the W and 
tt backgrounds are reduced by factors 3 — 6, while we maintain roughly 90 % of the 
BSM signal. 





no cut 


Minv = 

2 mr;() 
Al 


OGeV 


400 GeV 


10^ GeV 


-^vis, min 

[GeV] 
400 


500 


-^^cfT, min 

[GeV] 
400 


500 


W+ jets 


254.2 


62.53 


81.13 


48.73 


46.73 


70.29 


42.41 


75.89 


43.3 


tt+ jets 


151.4 


52.99 


64.67 


51.77 


50.14 


64.63 


46.3 


63.98 


40.12 


BSM 


31.05 


28.82 


29.82 


28.42 


27.80 


27.07 


24.43 


29.99 


28.72 



Table 9: Total cross sections a [pb] for BSM as well as VI^ + jets and tt + jets back- 
grounds, without and with several cuts on different inclusive quantities. First column: 
no cut; second to fifth column: values for y/Icut = ^inv + 500 GeV, with varying Minv 
values; last four columns: cuts on Mvis (Eq. (7)) and Mcs (Eq. (10)) respectively, where 
Minv = 2m^o. While the maximal suppression factor for the BSM signal is around 
1.11, the SM backgrounds are suppressed by factors 2 — 5. Equal results, however, can 
easily be obtained by a cut on Mvis or Mcs ■ 
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5.1 Comparison with other (transverse) variables 



The strength of V^min vs other (transverse) variables lies in the conjecture, given by 
Eqn. (6), about the direct correlation between its peak position (given the correct mass 
Minv) and the rise of the parton level production cross section. However, this correlation 
has so far not risen beyond the status of a conjecture. We therefore briefly discuss two 
other variables which might serve a similar purpose in background subtraction, namely 

Mvis = \/i^v'is-^?is (7) 

and 

$T = \PtI (8) 

where everything is defined at analysis object level. In [10], two more variables, namely 
Et and Ht = Et + $ti are studied^^; however, as we define 

= -'f^ (9) 

these two variables are only variations of and therefore not discussed here. Figure 
10 shows the subtracted (BSM+SM) — (SM) distributions of M^is and $j- respectively; 
especially the former looks quite promising. Indeed, a cut Minv > 500 GeV reduces 
the SM background by a factor 6 {W + jets) and 3 (tt+jets), cf. Table 9. We equally 
compare to the frequently used variable Meg [6, 37] 

vis 

which exhibits a similar power for background suppression as M^is, cf. Table 9. We 
therefore conclude that, for background suppression, both Mvis, Mcs as well as v^min 
work in a similar way; the advantage of the latter variable is the (conjectured) M^^^- 
dependent correlation between the its peak position and the hard (=parton level) 
process center-of-mass energy, if this can be proven to hold in all cases. Indeed, a 
similar conjecure of a linear correlation between the SUSY scale and the peak position 
of Mefj [6, 37] has recently been shown not to hold in all cases [38]; however, an 
equivalent systematic study of y/^mm is still lacking. 



6 Conclusion and outlook 

We investigated the variable \/§min for a complete BSM sample which includes all 
strong production as well as all possible decay chains. In our analysis, we include 
both soft as well as detector effects by including a complete parton shower as well as a 
generic detector simulation and reconstruction-level objects, which have been defined 
such that the parton-level variable can be recovered quite accurately. We investigate 
the variable V^min for a fully inclusive sample which sums over all possible final states 
of the hard scattering process, as well as for dominant exclusive final states. We see 
that, on parton level and for a correct input value of Mjnv, the y/^mm variable peaks 

^^Note that the definition of Ht differs in [13], where it basicahy is set to the variable Meff [6, 37]. 
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Figure 10: Difference between total and SM M^i^ (Eq. (7)) (left) and $t (Eq. (8)) 
(right) distributions for the total inclusive sample. No further cuts were applied. 



closely to the heavy particle production cross section which is in agreement with the 
conjecture made by the authors of [10]. In a comparison between parton level and anal- 
ysis level quantities, we see that in our sample the largest shift between these arises 
from the transition from tau-leptons which were used for the parton level quantity to 
the tau-jets and associated invisible neutrinos, which were used at analysis level. In 
order to asses this effect, we used the true tau-lepton four- vectors in the analysis level 
quantities. The effect is usually of the order of 100 GeV, and in the totally inclusive 
sample the shift completely disappears for the idealized case of perfectly reconstructed 
tau jets. We therefore conclude that for the parameter point considered here, even at 
analysis level, the parton level \/5min peak position can be sufficiently reconstructed^^. 
In case the correlation between the threshold of the parton level cross section and the 
peak of the V^min distribution could be proven rigorously, this would indeed provide a 
quite elegant and straightforward way to assess the scale of the new physics signal as 
a function of the total invisible mass of the process. 

Furthermore, we present the first study which investigates the use of V^min in order 
to suppress SM background for BSM searches. For this, we considered W+ jets as well 
as tt+jets background. We saw that these backgrounds could be sufficiently reduced 
by cuts on v^min, leading to suppression factors around 2 — 6, while we retained 90% of 
the BSM signal. This feature was independent of the input value of the total invisible 
mass Minv. However, we could achieve similar results by a cut on the total visible mass 
Mvis, which is a simpler variable which additionally does not require the input of Minv. 
A further comparison with M^s as a cut variable lead to similar results. We therefore 
conclude that, unless the conjecture about mass particle threshold and peak position 
of V^min can be rigorously proven, the latter does not exhibit significant advantages 
over other (transverse) variables. However, if the correlation between the threshold 
and "\/Imin could be rigorously proven, it would indeed provide a simple and elegant 
hold on the scale of new physics processes. A further investigation of this relation is in 
the line of future work. 

"'^^We point out that, for the scenario considered in this study, even when using analysis tau jets the 
maximal shift between parton level and analysis level peak positions was ^ 200 GeV, which effectively 
leads in an error ^ 100 GeV in the estimation of the initial heavy particles masses. The magnitude 
of this effect can of course differ depending on the BSM model and as well as specific model scenario 
point. 
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Appendix 

A SPSla spectrum 
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Table 10: Relevant masses for SPSla in GeV. u = (m,c), d = {d,s), I = {e,fi). 



B Minimization of \/^ 

As stated in Section 2, V^min denotes the minimal center of mass energy of an event 
with a measured visible four-vector P^^, = {E, 'Pti Pz) which is still in agreement with 
total energy momentum conservation as well as onshellness of all outgoing particles. 
We equally assume the event to be at rest in the transverse plane such that Eqn. (4) 
holds. The generic expression for s in this case is given by 

v..-^.) = - (11) 

where the index j goes over the invisible particles in the event with the respective 
energies 

7712 2 I 2 I 2 

Here and in the following, we omit the vector notation in px for simplification, but all 
transverse quantities should be read as pjr = {Pjx,Pjy), fr = {fxifv)) stc. 
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We use a Lagrange multiplier A to take the additional constraint for the vector sum of 
the transverse momenta 

YIpjT = fr (12) 

3 

into account. We therefore aim at minimizing 

£ = s 



iT — JTt 



ie we try to find the values of the invisible particles' three-momenta such that 

dC 



0. 



If we consider the case of riinv invisible particles in the event, we obtain 3 x ni, 
equations 



dC 

dpiT 

dC 

dpiz 



(13) 



Together with the constraint in Eqn. (12), we now have in total Srii^^ + 1 constraints 
for 3?T,inv + 1 unknowns {pix,Piy,Piz] A). From Eqns. (13), we immediately see that 



Ej _ PiT _ Piz 
Ej PjT Pjz 



Cij 



with Cij = ^ = const. Combining this with Eqn. (12), we obtain 

p,T ~ fr' 

where we defined M\^^ = rrij to be the sum over the masses of all invisible particles, 
c.f. Eqn. (2). We therefore have 

PiT = -T7 "Ij. 

^ M- 



We can now rewrite the second equation in (13) and obtain 

PiT J Ei \ PiT 



Solving this for piz leads to 



Pz'^nrii 



inv 



23 



We here have reproduced the solutions for pi^^, Pi,T given in [10] which minimize s. 
Inserting these into Eqn.(ll) then leads to VSmm given by Eqn. (1), which denotes the 
minimal hard scattering center of mass energy which is allowed by energy momentum 
conservation for a specific visible total four vector {E, P^, Pz) obtained from measure- 
ment. The only unknown quantity is Mjnv defined according to Eqn. (2), which has to 
be treated as an external input parameter for v'Imin- 



We want to comment that the transverse mass variable Mr [39, 40, 41, 42] has a 
functional form similar to V^min as given in Eqn. (1). For a system with visible and 
invisible total four- vectors P^g, P^^-, this variable is defined as 

with the transverse energies 

^T,vis = Mi + P|, El,,^ = (MW + fl. 

Here, Mvis, M[^^ denote the Lorentz-invariant masses of the total visible and invisible 
system respectively, 

VIS VIS' V mv/ mv' 

which vary on an event by event basis. Assuming Eqn. (4) to hold, it follows that 



We now see that the functional forms of Mt{MI^^) and "\/5min(Minv) as given in Eqn. 
(1) are identical, and differences between the variables only stem from the difference 
between M[^^ and Mjnv Indeed, for a correct guess of Mjnv, \f§ram and coincide if 
for the invisible particles in the event 

Componentwise, this equation is only fulfilled if we either have a complete set of mass- 
less particles which are all collinear with each other such that cos 6ij = 1 for all {i,j) 
pairs or for a complete set of massive particles which are all produced at rest. In 
general, however, 

ML > Mt:"'^ 



and therefore 



Mt > Vlmin ( 

on an event by event basis. 



^(truc)\ 



C Pythia 6.4 ISR/ FSR default setup 

All switch descriptions here are taken from [23]. We equally refer the reader to section 
10 of this reference for a more detailed discussion of the parton shower model and its 
implementation in Pythia. 
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MSTP(32) : (D = 8) definition in hard scattering for 2 — )■ 2 processes. For reso- 
nance production Q"^ is always cliosen to be s = m^, wliere itir is tlie mass of 
the resonance. 

The newer options 6-10 are specifically intended for processes with incoming 
virtual photons. These are ordered from a 'minimal' dependence on the vir- 
tualities to a 'maximal' one, based on reasonable kinematics considerations. 
The old default value MSTP(32) = 2 forms the starting point, with no de- 
pendence at all, and the new default is some intermediate choice. Notation 
is that and P2 are the virtualities of the two incoming particles, p± the 
transverse momentum of the scattering process, and and 7714 the masses 
of the two outgoing partons. For a direct photon, is the photon virtuality 
and X = 1. For a resolved photon, P^ still refers to the photon, rather than 
the unknown virtuality of the reacting parton in the photon, and x is the 
momentum fraction taken by this parton. 
= 2 : = (mis -|- m'i^)/2 = pl + {ml + m|)/2. 

= 8 : Q'^ = P± + {Pi + P2 +ml + m|)/2. ensure that the scale is always 
bigger than P^. 

MSTP(62) : (D = 3) level of coherence imposed on the space-like parton-shower evo- 
lution. 

= 3 : Q'^/p'± values and opening angles of emitted (on-mass-shell or time-like) 
partons are both strictly ordered, increasing towards the hard interaction. 
MSTP(63) : (D = 2) structure of associated time-like showers, i.e. showers initiated 
by emission off the incoming space-like partons in PYSSPA. 
= 2 : a shower may evolve, with maximum allowed time-like virtuality set by 
phase space or by PARP(71) times the value of the space-like parton 
created in the same vertex, whichever is the stronger constraint. 
MSTP(64) : (D = 2) choice of as and scale in space-like parton showers in PYSSPA. 

= 2 : first-order running as with argument PARP(64)/c^ =PARP(64)(1 - z)Q'^. 
MSTP(65) : (D = 1) treatment of soft-gluon emission in space-like parton-shower evo- 
lution in PYSSPA. 

= 1 : soft-gluon emission is resummed and included together with the hard 
radiation as an effective z shift. 
MSTP(66) : (D = 5) choice of lower cut-off for initial-state QCD radiation in VMD or 
anomalous photoproduction events, and matching to primordial k±. 

= 1 : for anomalous photons, the lower cut-off is the larger of PARP(62)^ 
and VINT(283) or VINT(284), where the latter is the virtuality scale for 
the 7 — )■ qq vertex on the appropriate side of the event. The VINT values 
are selected logarithmically even between PARP(15)^ and the scale of 
the parton distributions of the hard process. 

= 4 : a stronger damping at large k±, like dk\/{k'j_ + with ko < kj_ < 

P±mm{W^). Apart from this, it works like = 1. 

= 5 : a k± generated as in = 4 is added vectorially with a standard Gaussian 
k± generated like for VMD states. Ensures that GVMD has typical k±^s 
above those of VMD, in spite of the large primordial k±^s implied by 
hadronic physics. (Probably attributable to a lack of soft QCD radiation 
in parton showers.) 
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MSTP (67) : (D = 2) possibility to introduce colour coherence effects in the first branch- 
ing of the backwards evolution of an initial-state shower in PYSSPA; mainly of 
relevance for QCD parton-parton scattering processes. 
= 2 : restrict the polar angle of a branching to be smaller than the scattering 

angle of the relevant colour flow. 
Note 1: azimuthal anisotropics have not yet been included. 

Note 2: for subsequent branchings, MSTP (62) = 3 is used to restrict the (polar) 
angular range of branchings. 
MSTP (68) : (D = 3) choice of maximum virtuality scale and matrix-element matching 
scheme for initial-state radiation. To this end, the basic scattering processes 
are classified as belonging to one or several of the following categories (hard- 
coded for each process): 
= : maximum shower virtuality is the same as the choice for the parton 
distributions, see MSTP (32). (Except that the multiplicative extra factor 
PARP(34) is absent and instead PARP(67) can be used for this purpose.) 
No matrix-element correction. 
= 3 : as = 0, but ME corrections are applied where available. 
MSTP (69) : (D = 0) possibility to change scale for parton distributions from the 
MSTP (32) choice, especially for e+e". 
= : use MSTP (32) scale. 
MSTP (72) : (D = 1) maximum scale for radiation off FSR dipoles stretched between 
ISR partons in the new p^-ordered evolution in PYPTIS. 
= 1 : the p_Lmax scale of FSR is set as the p_L production scale of the respective 
radiating parton. Dipoles stretched to remnants do not radiate. 

The additional switches/ variables appearing above are given by 

PARP(15) : (D = 0.5 GeV) lower cut-off po used to define minimum transverse mo- 
mentum in branchings 7 — )■ qq in the anomalous event class of 7p interactions, 
i.e. sets the dividing line between the VMD and GVMD event classes. 

PARP(62) : (D = 1. GeV) effective cut-off Q or /c_l value (see MSTP(64)), below which 
space-like parton showers are not evolved. Primarily intended for QCD show- 
ers in incoming hadrons, but also applied to q — ?■ q7 branchings. 

PARP(64) : (D = 1.) in space-like parton-shower evolution the squared transverse 
momentum evolution scale is multiplied by PARP(64) for use as a scale in 
and parton distributions when MSTP (64) = 2. 

PARP(67) : (D = 4.) the scale of the hard scattering (see MSTP (32)) is multiplied 
by PARP(67) to define the maximum parton virtuality allowed in Q^-ordered 
space-like showers. This does not apply to s-channel resonances, where the 
m aximum virtuality is set by m^. It does apply to all user-defined pro- 
cesses, however. 

PARP(71) : (D = 4.) the scale of the hard scattering (see MSTP (32)) is multiplied 
by PARP(71) to define the maximum parton virtuality allowed in time-like 
showers. This does not apply to s-channel resonances, where the maximum 
virtuality is set by . Like for PARP(67) this number is uncertain. 

VINT(283), VINT(284) : virtuality scale at which a GVMD/anomalous photon on 
the beam or target side of the event is being resolved. More precisely, it gives 
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the k\ of the 7 — qq vertex. For elastic and diffractive scatterings, is 
stored, where m is the mass of the state being diffracted. For clarity, we point 
out that elastic and diffractive events are characterized by the mass of the 
diffractive states but without any primordial fc_L, while jet production involves 
a primordial kj_ but no mass selection. Both are thus not used at the same 
time, but for GVMD / anomalous photons, the standard (though approximate) 
identification k'j_ = ensures agreement between the two applications. 

VDM/ GVDM are acronyms for vector meson dominated/ generalized vector meson 
dominated events in photo production respectively (cf section 7.7.2 of [23]). 
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